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Abstract. For the class of polynomial skew products of C , we describe a 
rigorous computer algorithm which, for a given map /, will (1) build a model 
of the dynamics of / on its chain recurrent set, and (2) attempt to determine 
whether / is Axiom A. Further, we discuss how we used our implementation 
of this algorithm to establish Axiom A for several explicit cases. 



1. Introduction 

Our main interest here is to develop and use rigorous computer investigations 
to study the dynamics of polynomial skew products of C 2 ; i.e., maps of the form 
/(z, w) = (p(z), q(z, w)), where p and q are polynomials of the same degree d > 2. 

The skew products we are most interested in studying are those maps which are 
Axiom A. Such maps have the "simplest" chaotic dynamics, and stability under 
small perturbation, thus are amenable to computer investigation. In one complex 
dimension, a polynomial map is called hyperbolic if it is uniformly expanding on 
some neighborhood of its Julia set, with respect to some riemannian metric. In 
C n , a map is hyperbolic on an invariant set A if there exists a continuous splitting 
of the tangent bundle over A into two subspaces (of any dimension zero through 
n), with one subspace uniformly expanded by the map, and the other uniformly 
contracted. A map is Axiom A if / is hyperbolic on the nonwandering set, f2, and 
if periodic points are dense in f2 and f2 is compact. 

Our motivation is to understand what kind of dynamics can occur for Axiom 
A polynomial skew products. In this paper, we construct a class of maps with 
interesting dynamics, and prove using rigorous computer techniques that sample 
maps from this class are Axiom A. This leads us to conjecture that all (or nearly 
all) maps in this class are Axiom A. 

To develop our computer techniques, we use the work of |l(Jj and as a 
foundation. In |l(Jj . we described a rigorous algorithm (and its implementation) 
for constructing a neighborhood 23 of the chain recurrent set, 31, and a graph T 
modelling the dynamics of / on 23. In j 1 1 j we developed a numerical method for 
proving hyperbolicity of a polynomial map p of C, which consists of an attempt to 
construct a metric in which p is expanding on 23, by some uniform factor L > 1. 

Here, we adapt the one-dimensional algorithm to the setting of skew products. 
Our efforts are aided by the fact that skew products are a natural generalization of 
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polynomials in one dimension, since a skew product / maps the vertical line {z} x C 
to the vertical line {p(z)} x C, and restricted to a vertical line, / is the polynomial 
map w i— > q z (w) = q(z, w). 

We applied a similar approach in where we developed a test for verifying 
hypcrbolicity of polynomial diffcomorphisms of C 2 . This algorithm is different 
from that for skew products, since the dynamics of polynomial diffcomorphisms is 
of saddle type, i.e., with one direction expanding and the other contracting, whereas 
for skew products, it happens that to establish Axiom A we need only check an 
expanding condition. 

We have implemented all of our algorithms into a computer program called 
Hypatia, which began with the algorithms of ^Oh then was enhanced to include 
algorithms from (and [H]), and now with the present work, can be used to 
prove Axiom A for specific polynomial skew products. 

Finally, we describe the organization of this paper. Polynomial skew products of 
C 2 have been studied by Jonsson f|12ll5|') and Heinemann (00111]). In Section |2 
we provide necessary background on skew products, hypcrbolicity, and invariant 
sets of interest. In Section 03 we describe the dynamics of several classes of skew 
products, including our new class discussed above. In Section 0] we outline our 
rigorous computer algorithms, implemented in Hypatia, for attempting to verify 
Axiom A for skew products. In Section [5] we describe how we used Hypatia to 
prove Axiom A for several specific skew products, including maps of the same 
type as Jonsson's and Hcincmann's examples, in addition to our new example. In 
Appendix EI we give a brief overview of the technique that we used to control 
round-off error in our computations, called interval arithmetic. 

Acknowledgements. I would like to thank John Smillie for suggesting this project, 
Mattias Jonsson and Manfred Dcnkcr for advice on getting started, and Eric Bed- 
ford for guidance throughout my investigations. I would also like to thank Adricn 
Douady and Mikhail Lyubich for asking enlightening questions when I presented 
this work at the Snowbird meeting in June 2004, and Amie Wilkinson for mention- 
ing the relationship between polynomial skew products and partial hypcrbolicity. 

2. Background 

2.1. Chain recurrence and hyperbolicity. The chain recurrent set, 3i, the Julia 
set, J, and the non-wandering set, Q, are all attempts at locating the points with 
dynamically interesting behavior. The Julia set, J, can be defined as the topological 
boundary of the set, K, of points in C 2 with bounded orbits under /. Both J and 
K are compact. Slices of the Julia set can be easily sketched by computer, but 3? 
is the set most amenable to rigorous computer investigation. 3i can also be easily 
decomposed into components which do not interact with one another. Since 'R D J, 
we can learn about J by studying / on 3?. 

An e-chain of length n > 1 from y to z is a sequence of points {y = x\, . . . ,x n = 
z) such that \f(xk) — %k+i\ < e for 1 < i < n — 1. A point y belongs to the e- 
chain recurrent set, 3i £ , of a function / if there is an e-chain from y to y. The 
chain recurrent set is 3? = n e >o3l £ . A point z is in the forward chain limit set of 
a point y, 3i(y), if for all e > 0, for all n > 1, there is an e-chain from y to z of 
length greater than n. Put an equivalence relation on 5i by: y ~ z if y G 3i(z) and 
z e 3£(y). Equivalence classes are called chain transitive components. Define Oi e (y) 
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and e-chain transitive components analogously. ft is closed and invariant, and if 
e < e', then ft C ft e C ft e *. 

Thus chain recurrence is quite natural to rigorously study using a computer. 
Next, we recall the precise definitions of hypcrbolicity and Axiom A. 

Definition 2.1. Let g be a C 1 diffcomorphism or endomorphism of a compact 
manifold M, and let A be a closed, ^-invariant set. Say A is hyperbolic for g if there 
is a splitting of the tangent bundle T X M = E s x © E% (one subspace may be trivial), 
for each x in A, which varies continuously with a; in A, a constant A > 1, and a 
riemannian norm ||-|| such that: 

(1) g preserves the splitting, i.e., D x g{E s x ) = Eg X , and D x g{E x ) = E^ x , and 

(2) Dg expands (contracts) E U {E S ) uniformly, i.e., ||£) x g(w)|| > A||w|| for 
w G El, and ||D x ff(v)|| < (1/A) ||v|| for v G E%. 

This definition is independent of choice of norm. 

Definition 2.2. g is Axiom A if for the nonwandering set f2, we have: (1) fl is 
compact, (2) periodic points are dense in £1, and (3) f2 is hyperbolic for g. 

If Pcr(g) is the set of periodic points of g, then Per(g) C fi C ft. Jonsson shows 
that for Axiom A polynomial skew products, Per(/) = fi = ft (see Theorem 12.311 . 

2.2. Polynomial skew products of C 2 . In this subsection, we summarize some 
of the notation and results of to give needed background on skew products. 

Since the dynamics of / in the z-coordinate is given by p, it will be useful to 
employ the notation (K p and) J p for the one-dimensional (filled) Julia set of p, and 
G p (z) for the Green function in C of p, where K p = {G p = 0}. 

Let M denote the Mandelbrot set for quadratic polynomials, i.e., M is the set of 
all c such that the critical orbit is bounded under p c (x) = x 2 + c (equivalently, M 
is the locus of c with connected J c ). 

Global Dynamics. For polynomial skew products, the usual rate of escape 
Green function, defined for x G C 2 by G(x) = lim n ^oo -t log + |/(x)|, is continuous, 
plurisubharmonic, nonnegativc, and satisfies Go / = dG and K = {G = 0}. One 
can also define a positive closed current T = ^dd c G and an ergodic invariant 
measure, \i = T A T, of maximal entropy logd 2 . J2 := supp(/i) is also the closure 
of the set of repelling periodic points. 

Vertical dynamics. Since / preserves the vertical lines {z} x C, it is useful to 
consider the dynamics of / on this family of lines. Let z n — p n (z), q z (w) — q(z, w), 
and Q"(u>) = q Zn _ 1 o- • - oq z , so that f n (z, w) = (z n , Q™(w)). Let G z (w) = G(z,w)- 
G p (z). Then G z is nonnegative, continuous, subharmonic, and is asymptotic to 
log|w| — G p (z) as w — * 00. Naturally, define K z = {G z = 0}, and J z = dK z . 
Then K z and J z are compact, and if z G K p , then w G K z if and only if |Q™(w)| is 
bounded. Further, G Zl o q z = dG Zl which implies q z {K z ) = K Zl and q z {J z ) = J Zl . 
Define C z as the critical set in {z} x C. For example, if q(z, w) = w 2 + b + ez, then 
C z = {0} for every z. For all z G C, J z is connected if and only if C Zn C K Zn . 

However, not every phenomena of one-dimensional dynamics carries over to ver- 
tical dynamics. For example, unlike in one dimension, J z may have finitely many 
(but greater than one) connected components, even for d = 2 (see |12j . remark 2.5). 

Vertical Expansion. Let Z C K p be compact with pZ C Z , for example 
Z = J p or Z = A p , the set of attracting periodic orbits. Let Jz = ] ^ z <ez{z} x J z . 
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Jonsson shows Jj = J 2 - Call / vertically expanding over Z if there exist c > and 
A > 1 such that \DQ n z (w)\ > c\ n , for all z G Z, w G J z , and n > 1. 

Let Cz = ^zez{z}xC z and Dz = U„>i/™Cz- Further, / is vertically expanding 
over Z if and only if Dz C\ Jz = 0- 

As in one dimension, z 1— ► A", is upper semicontinuous and z 1— > J z is lower 
semicontinuous, in the Hausdorff metric. Further, if / is vertically expanding over 
Z , then z t— ► J z is continuous for all zeZ; and if in addition, J z is connected for 
all 2 G Z, then z 1— > A z is continuous for all z G Z. However, Jonsson provides 
examples showing vertical expansion over J p U A p neither implies that z 1— > J z is 
continuous on all of C, nor that z t— > K z is continuous for z G J p . We will study 
one example of the latter phenomenon in Section |SJ 

Axiom A polynomial skew products. In this paper we rely significantly on 
Jonsson's vertical expansion criteria for Axiom A: 

Theorem 2.3 (^2L Theorem 8.2). A polynomial skew product f is Axiom A on 
C 2 if and only if 

(1) p is uniformly expanding on J p , 

(2) / is vertically expanding over J p , and 

(3) / is vertically expanding over A p . 

Moreover, if f is Axiom A, then 31 = Clf = Per(f). 

Note Jonsson's criticality condition for vertical expansion over Z implies 

Corollary 2.4 (|12|. Corollary 8.3). A polynomial skew product f is Axiom A on 
C 2 if and only if (1) D p n J p = 0, (2) D Jp H J 2 = 0, and (3) D Ap n J Ap = 0- 

Jonsson also provides a structural stability result for Axiom A skew products 
(\V2\. Theorem A. 6 and Proposition A. 7). A consequence is that it makes sense 
to refer to a connected component of the subset of Axiom A mappings in a given 
parameter space as a hyperbolic component. 

3. Examples of Axiom A polynomial skew products 

Products. A simple product, f(z, w) = (p(z), q(w)) is Axiom A if p and q are 
hyperbolic. There are (at most) four chain transitive components: 51 = (J p x J q ) U 
(A p x J q ) U (J p x A q ) U (Ap x A q ); the first is J2, the codimension zero set on which 
/ is uniformly expanding, the middle two are codimension one saddle sets, and the 
last set is the attracting periodic points of /. For example, for (z,w) 1— > (z 2 ,w 2 ), 
we have J p x J q = S 1 x S* 1 is a torus, A p x J q = {0} x S 1 and J p x A q = S 1 x {0} 
are circles, and A p x A q is the origin. 

(Small) Perturbations of products. By Jonsson's structural stability re- 
sults, any f(z, w) which is a sufficiently small perturbation of an Axiom A product 
is also Axiom A, with basic sets (i.e., chain transitive components) topologically 
corresponding to those of the product. 

Jonsson's nonproduct example. In [5], Diller and Jonsson describe a class of 
maps which are Axiom A, but not conjugate to any Axiom A product (again using 
Corollary I2.4[l . and consider the example: fo(z,w) = (z 2 — 9900, w 2 + (z + 99)/6). 

This type of example starts with p(z) — z 2 — R with R > 2, so J p is a real Cantor 
set. Let a < < j3 denote the two (necessarily repelling) fixed points of p. Let 
D x := [-/?, -17], where 77 = y/R-fi, and let D 2 = -D x . Then J p C D x U D 2 . For 
fo, we have R = 9900, (3 = 100, a = -99, and n = V^OO « 98.99. 
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The skew product /o is of the form (z 2 — R, w 2 + (z + f3)/S), where S > 1. 
Thus for z G £>i, i.e., close to — /?, we have q z {w) ~ u> 2 . One can easily check that 
a + [3 = 1, so if S is sufficiently large, then q a (w) — w 2 + (a + (3)/S is conjugate 
to w 2 , and so since a is a fixed point, J Q is a quasicircle. On the other hand, for 
z G £>2, i-e., z close to f3, we have q z (w) « u> 2 + 2(3/S, and since /3 > 2, if 5 is 
not too large, then q z (w) m w 2 + c for some c > outside the Mandelbrot set, 
i.e., with escaping critical point. Hence Jp is a Cantor set. Using the fact that 
q(J z ) = Jp( z ), we quickly see that the fibers J z for z G J p contain a rich mix of 
Cantor sets and collections of circles. The nonwandering set consists of J 2 (the 
closure of the repelling periodic points) and the saddle fixed point (a,0). If R is 
sufficiently large, then such a map is Axiom A. 

Dillcr and Jonsson describe a method for generating an Axiom A example of this 
type for every degree d > 2 (in fact, for every k with 1 < k < d, they can construct 
an example with topological entropy log k) . J p is always contained in the union of 
two disjoint sets, E and E\ with q z (w) w d for z G E and q z (w) = w d + r(z), 
such that r(z) is large for z G E' . but r(z) is of degree less than d. 

Generating new examples. The first proposition below is an attempt to 
generalize the example of Diller and Jonsson, by replacing the map w 2 , and the 
circles in J that arise from this map, with the map q c {w) = w 2 +c for any hyperbolic 
c G M, and thus replacing circles with any connected hyperbolic Julia set J c . 

Proposition 3.1. For any hyperbolic polynomial x i— > x 2 + c, with c G M, there 
exist a G C, R > 6, and S > such that for 



(1) J p C D\ U D 2 , for the intervals D\ = [—[3, —rj\ and D 2 = —D\, where 
(3 > > a are the fixed points of p, and rj = y'R — f3; 

(2) for all z G D\ , q z is in the same hyperbolic component of M as c (hence 
q z (w) is topologically conjugate to w 2 + c); and 

(3) for all z G D2, q z is outside o/M. 

Proof. First, we quickly check (1). The fixed points a, (3 are given by: a = (1 — 
\/l + 4R)/2,[3= 1 + \/l + 4R)/2. If R > 2, then straightforward calculation gives 
R > (3 > 2, hence 77 > 0, and p _1 ([-A =DiU L> 2 , with p _1 (/3) = ±/3 and 
p _1 (— /?) = ±r/ (and p _1 (a) = ±a), and p a homcomorphism on each of D\ and 
D 2 . The map p(z) = z 2 — R is conjugate to the one-sided two shift, with J p a real 
Cantor set contained in D\ U D 2 - 

Let a = (77 + /3)/2. Now, given a, we find an S* establishing (2), for any R > 6. 
Note over all the z G Z?i, the interval of g z 's is exactly centered around w 2 + c, with 
the extreme q z s differing from w 2 +c by length(_Di)/(25) = ((3 — rj)/(2S). We just 
need to require S to be large enough to keep q z in the same hyperbolic component 
as w 2 + c. Let a denote the horizontal distance between c and the boundary of the 
hyperbolic component of M containing c. So we need ((3—r])/(2S) < a. Since R > 6, 
[3-r, = l3- < 3 - V3. Hence if S > 2<r/(3 - V3), then {(3 - rj)/(2S) < a. 

Finally, we show that given a and S, we can pick R large enough to satisfy (3). 
For all z G D 2 = [?/, /3], we need <? z (iu) w w 2 + large, i.e., we need g z outside of 
M. The smallest q z is q v (w) = w 2 + c + (rj + a)/S = w 2 + c+ (3r/ + (3)/(2S). So 
we just need Re(c) + (3r/ + (3)/(2S) > 2 to guarantee we're outside of M. Since 




z 2 - R, w 2 +c + 




we have: 
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ceMwe know Re(c) > —2, so we just need (3?? + (3)/(2S) > 4. But S is fixed, so 
all we must do is make (3?y + (3) large, which we can do by taking R large enough. 
Specifically, let R be large enough that ?/ > 25. We can do this since r} <= ^/R — [3 
and (3 ~ VR, so 77 tends to infinity as R does. Since (3 > 17, this yields the required 
(3tj + (3) (25) > 4?7/(2S*) = 2^/S > 2(2S)/S = 4. □ 

Remark 1. We conjecture that all, or at least a large subset of, maps satisfying 
Proposition 13 . II are Axiom A. 

Unfortunately, we do not know a general method for verifying the above remark. 
We prove Axiom A for a specific example using the program Hypatia in Section 

One would ideally generalize this class of examples further, by replacing the 
Cantor sets above with an arbitrary hyperbolic polynomial, i.e., construct a skew 
product of the form (z 2 — R,w 2 + l(z)), with l(z) some linear or quadratic map, but 
instead of (3) above require q z for z G D2 to be in the same hyperbolic component 
of M as some other hyperbolic map x 2 + C2 ■ 

However, this does not work in general for any two hyperbolic maps x 2 + c\, x 2 + 
C2- For, we would need l(z) to map [—/?,/?] into the parameter space of c's, so that 
l(Di) is contained in the hyperbolic component of M containing c,. But R > 2 is 
necessary, and so the length of Di is not too small, but the hyperbolic components 
containing c, may be small compared to the distance between ci and C2. Hence 
there may be no linear or quadratic map l(z) which will work. 

Hence the best we can do is consider ci,C2 G M hyperbolic, ideally near the 
center of their hyperbolic components, then choose some R > 6 and consider 

(1) f(z, w) = (p(z), q z (w)) =(z 2 -R, w 2 + |(l - J) + ! (l + -)) , 

where a = (?] + /3)/2 is the center of D 2 as in Proposition l3.il Then q z (w) is of the 
form w 2 +l(z), and l(z) is the linear map satisfying ?(— a) = C\ and 1(a) = c 2 . Then 
we must check whether l{—/3) and ? ( — 97) are in the same hyperbolic component as 
C\, and whether l(rj) and l((3) are in the same hyperbolic component as 02- If this 
holds, the i(-Dj) is in the same hyperbolic component as Cj, and we have the desired 
conclusion. In this case, we conjecture the map is Axiom A. 

Fortunately, for the several largest hyperbolic components of M, the above con- 
struction is successful. For example: 

Example 3.2 (Circles and Basillicas). For f(z,w) = (p(z),q z (w)) = (z 2 — 90, u> 2 — 
z/20 - 1/2), it is easy to check that: (1) J p C £>i U D 2 , with D 1 = [-10, ~V80}; 

(2) for all z E D%, q z (w) is in the same hyperbolic component as w 2 ; and (3) for all 
z G -D2, Qz(w) is in the same hyperbolic component as w 2 — 1 (the Basillica). 

Example 3.3 (Basillicas and Rabbits). For f(z, w) = (p(z),q z (w)) = (z 2 — 90, w 2 + 
z(0.046 + 0.04z) - .56 + .37i), it is easy to check that: (1) J p C D x U D 2 , with 
D\ = [—10, — VSO]; (2) for all z G D\, q z (w) is in the same hyperbolic component 
as w 2 — 1 (the Basillica); and (3) for all z G D 2 , q z (w) is in the same hyperbolic 
component as w 2 — 0.12 + 0.75i (the Rabbit). 

4. Hypatia FOR POLYNOMIAL SKEW PRODUCTS 

In this section wc describe the algorithms for studying polynomial skew products 
that we have implemented in Hypatia. There are two main steps in Hypatia: 
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(1) building a neighborhood 23 of some 51$ and a graph modelling the 1-step 
dynamics of / on 23, and 

(2) testing whether a skew product / is Axiom A, using Jonsson's vertical 
expansion criteria, by attempting to build a metric for which / is vertically 
expanding over neighborhoods of J p and A p . 

4.1. Building a model of the chain recurrent set. In ^H] we describe an algo- 
rithm we call the box chain construction, for building a directed graph T modelling 
a given map g of C' 1 on a neighborhood 23 of 51. A similar approach in different 
settings can be found in [3 El El El H61 ITT| . The basic output of this construction 
is a graph called a box chain recurrent model, which we define below. 

Definition 4.1. Let A be an invariant set of a map g: C™ — > C n . Let T = (V, £) 
be a directed graph, with vertex set V = {-Bfc}£L 1; a finite collection of closed boxes 
in C n , having disjoint interiors, and such that the union of the boxes 23 = U% =1 B k 
contains A. Suppose there is a S > such that T contains an edge from B k to Bj 
if the image g(Bk) intersects a ^-neighborhood of Bj, i.e., 

8=>{(*,i): fl (B fc )n>T(B Jj *) 9 fe0}. 

Further, assume Y is partitioned by edge-connected components T' which are strongly 
connected, i.e., for each pair of vertices Bk, Bj, there is a path in V from B k to Bj, 
and vice-versa. Then we say Y is a box chain recurrent model of g on A, and each 
r' is a box chain transitive component. 

A box chain recurrent model Y is an approximation to the dynamics of g on 51, 
and the connected components Y' of Y are approximations to the chain transitive 
components. 

The basic box chain construction for a polynomial map p of C can be summarized 
as follows: 

(1) Compute R such that for some 5 > 0, % C 23o := [-R, R] 2 - 

(2) Subdivide 23 into a grid of boxes {Bk}. 

(3) Build a graph T with vertices {Bk] and edges {(k,j) : p(Bk) fl N(Bj, S) ^ 
0}. 

(4) Find the maximal subgraph T of T which consists precisely of edges and 
vertices lying in cycles. Then Y is partitioned by its edge-connected com- 
ponents: r = (r 1 u r 2 u . . . u r^). 

(5) If desired, refine Y by subdividing the boxes of Y and repeating (3) and (4). 
Then T is a box chain recurrent model of p on 51, and the Y J are the box chain 
transitive components. 

Now we outline the process we follow for skew products. 

(1) Use the one dimensional box chain construction of ^0] to build a box chain 
recurrent model Y for p in the z-plane. Iterate the construction until the boxes are 
small enough that (a) Y separates the box chain transitive component containing J p , 
call it Y°j , from the one containing A p (the attracting periodic orbit, if there is one), 
called and (b) hyperbolicity of p can be established for r ° Jp (using the techniques 
of |11|). Construction begins by calculating an i?i such that 5i p c Ri) 2 . The 
boxes of r, {-Bfc}, are a subset of the boxes in a 2™ x 2™ grid on [— Ri] 2 ■ 

(2) Build a model of / in only the box fibers over Y° A and Y°j . 

That is, first compute an R 2 so that 5ls(f) C [— Ri,R{\ 2 x [— R2,R 2 ] 2 (this 
calculation is similar to computing Ri, which is explained in |10|). Then choose 
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some m > and construct a grid of 2 m x 2 m boxes {-B™} 2 =1 in the square [— R2, R2] 2 
in the w- plane. Form boxes in C 2 as the product of boxes in F^ or F j times boxes 
in the new ui-grid, i.e., for each box B%. in T A or T° , we have 2 m x 2 m boxes in 
C 2 of the form B kJ =B%x Bf . Thus we get a set of boxes in C 2 , {Bi}. 

Now build a transition graph for /, T, whose vertices are the boxes Bi and 
with an edge B k -> B 3 if f(B k ) n N(B 3 ,S) ^ 0. Finally, find the subgraph T of T 
consisting precisely of the vertices and edges of T which lie in cycles, and decompose 
r into its edge-connected components: T = (T 1 U T 2 U . . . U T N ). Each T J is a box 
chain transitive component. If the boxes in C 2 are sufficiently small, then the 
box chain transitive components separate the chain transitive components of /; for 
instance, any box chain transitive components over should be distinct from 
those over F^ . There may be more than one box chain transitive component over 
an invariant set in the base, for example if q(z, w) = w 2 , then there is a component 
containing S 1 = J z in the w-plane, and (if boxes are small enough) a separate 
component containing the attracting fixed point of q z , the iw-plane's origin. 

Recall Jz = l J z ^z{z} X Jz- Let denote the box chain transitive component 
over which contains Ja p ] similarly let r 2 ^ be the box chain transitive compo- 
nent over Tj containing Jj p . Heuristically, we can quickly identify the box chain 
transitive component over a set Z which contains Jz , because it will be the largest 
by far (i.e., have the most vertices and edges). 

4.2. Establishing Axiom A: box vertical expansion. After the modified box 
chain construction of the previous section has produced box chain transitive com- 
ponents (containing Ja p ) and F 2 (containing Jj p ), it is time to test whether 
/ is Axiom A. We use the vertical expansion criteria of Theorem l2.3l i.e., we want 
(1) p uniformly expanding on J p , (2) / vertically expanding over J pi and (3) / ver- 
tically expanding over A p . Note (1) was checked in the construction of the previous 
subsection, so in this subsection we deal with (2) and (3). 

Recall that / is vertically expanding over Z if there are c > 0, A > 1 such 
that \DQ™(w)\ > cA™, for all z £ Z, w G J z , and n > 1, where let z n = p n (z), 
q z (w) = q(z,w), and Q z (w) = q Zn _ t o • • - oq z , so that f n (z,w) = (z n , Q™(w)). Thus 
vertical expansion is expansion in one complex dimension only, specifically, in the 
w-coordinate. This allows us to easily adapt our algorithm from jllj for verifying 
hyperbolicity of polynomial maps of p to testing vertical expansion. 

In one dimension, hyperbolicity of p reduces to uniform expansion of tangent 
vectors over a neighborhood of J, in some riemannian metric. Hypatia tests hyper- 
bolicity by trying to build a piecewise constant metric in which p is expanding: 

Definition 4.2. Let p be a polynomial map of C. Let V be the box chain transitive 
component containing J p , with box vertices {Bk} k=1 . Call p box-expansive on T' if 
there exists L > 1 and positive constants {<pk}k=i sucn that for every edge (k,j) 
in r', we have <fij \D z p(v)\ > L<fik \v\, for all z £ Bk, and all v £ T Z C (or simply 
<f> 3 \p'(z)\ > Lcp k ). 

The constants {4>k}k=i define what we call a box metric on H>' = U^ =1 Bk, i.e., 
simply a piecewise constant multiple of the euclidean metric on 3', with a different 
constant for each box. In [llj . we establish the following theorem, showing that 
box expansion implies the standard notion of hyperbolicity: 
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Theorem 4.3. Suppose p is box-expansive by L> 1 on 23' = U% =1 B k - Let Ai = 
23'np _1 (23'). Then there exists an r\ > 1 and a continuous norm \-\ on Ta x C which 
p expands by r), i.e., \D z p(v)\ > 77 |u| , /or aZ/ 2 G Ai, and all v £ T Z C 

The notion of box expansion is easily adapted to vertical expansion. 

Definition 4.4. Let f(z,w) = (p{z) , q(z , w)) be a polynomial skew product. Let 
Z C K p be compact with pZ C Let L' be the box chain transitive component 
containing Jz = U z ^z{z} x J z , with box vertices {Bk}^ =1 . Call / box vertically 
expansive on L' if there exists L > 1 and positive constants {<^fc}fc=i such that for 
every edge (k,j) in L', we have </>j \D w q z (v)\ > Lfa |u|, for all (2, i/j) € Bfc, and all 
ner m C (or simply <j)j \q' z (w)\ > L<f) k ). 

Theorem 14.31 immediately yields as a corollary its analog for vertical expansion. 
We use that the definition of vertical expansion over Z given by Jonsson is equiv- 
alent to the existence of some constant 77 > 1 and some riemannian norm |-| in 
the tangent bundle to the w-coordinate over Z, such that in this norm, DQ Z ex- 
pands tangent vectors by rj. That is. if Jz C 23' = Ug^r'^t, then box vertical 
expansion for L' implies vertical expansion of / over a neighborhood of Z, in some 
(continuous) riemannian norm, on the u>-coordinate tangent bundle for all [z, w) in 
a neighborhood Ai, with Jz C Ai C 53'. 

We describe an algorithm for attempting to establish box expansion in ^J. To 
test box vertical expansion, we can use this algorithm with only very small changes. 
We perform the one-dimensional algorithm twice, once for L' = T\ (containing 
Ja p ) and then again for V = (containing Jj p ). For each L', we use the one- 
dimensional recursive algorithm to attempt to build a set of metric constants {4>k\ 
for each box Bk of F, such that for each edge (k,j) in L', the inequality 



W3 ~ mf \q' z (w)\ 

is satisfied. Note if any of the derivatives q' z (w) is zero, then L' contains a crit- 
ical point (of q z ), so box vertical expansion necessarily fails and the algorithm is 
terminated. At this point, the user may re-try with smaller boxes, to attempt to 
separate this critical point from Jz- If all such derivatives are positive, the recur- 
sive algorithm to define the (f>k may still fail. This will happen if the map is not 
vertically expanding over Z ', or it may mean the boxes or L are simply too large. 
Thus the user may re-try with smaller L, or smaller boxes. On the otheer hand, if 
the algorithm is successfully completed, producing a set of metric constants {</>&} 
for r', then Theorem 14.31 implies that / is vertically expanding over J p and A p , i.e., 
(2) and (3) of Theorem 12. 31 are satisfied, and combined with (1), which we checked 
in Section l4~Tl this yields that / is Axiom A. 



5. Results from Hypatia 

In this section, we describe the results of applying the algorithms of Section 0] 
to some specific examples of maps from the classes discussed in Sectional 

The computations described in this section were run via a C ++ program in a 
Unix environment, with 2GB of RAM. When computations became overwhelming, 
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memory usage was the limiting factor. Rigor was maintained in the computations 
using the Interval Arithmetic routines in the PROFIL/BIAS package, available at 
|18| . See Appendix lAl for a brief introduction to this technique. 

The examples we examine are all quadratic skew products. Thus p and q z each 
have only one critical point, at the origin, and p has at most one attracting cycle. 

Before proceeding with the examples, we summarize the data given. For each 
example, the box chain transitive components were formed by boxes from a (2 ra ) 2 x 
(2 m ) 2 grid on 23 = [— -Ri] 2 x [— R2,R-2} 2 - For each graph component for which 
vertical expansion was tested, we give the size of the graph, (V;E), with V the 
number of vertices and E the number of edges. If expansion was established, we 
state the expansion constant, L, the range of metric constants 4>k constructed, as 
an interval [0 m i n , </> ma x] , and the average metric constant </> a vg- 

Example 5.1 (Circle cross circle perturbation). The map f(z,w) = (z 2 ,w 2 + 
iu/lO + z/100) appears to be a stable perturbation of the product (z,w) i— ► (z 2 ,w 2 ). 
This map is in Heinemann's family of cannelloni maps, since J p is connected and 
the fiber Julia sets are quasicircles. 

We proved with Hypatia that it is Axiom A, using boxes from a (2 5 ) 2 x (2 4 ) 2 
grid on 23 = [—1.1, l.l] 2 x [—1.21, 1.21] 2 . Since z i— > z 2 has an attracting fixed point 
at the origin, there are two box chain transitive components to test expansion over: 
one for Jj p and one for Ja v = {0} x J{o}- The component containing Jj p , of 
size (V;E) = (39,200; 1,953,900), is box vertically expansive for L = 1.29105, with 
metric constants in [0.1674, 1], with average 0.290. The component containing J a , 
of size (V;E) = (608; 4,416), is box vertically expansive for L = 1.29105, with 
metric constants in [0.11116, 1], with average 0.168. 

Example 5.2 (Cantor cross circle perturbation). For the map f(z,w) = (z 2 + 
2, w 2 + z/10), we have J p is a Cantor set in the base, and q z is a small perturbation 
of ro n to 2 . Heinemann calls this type of perturbation a Cantor skew, and shows 
for sufficiently small perturbation, the fiber Julia sets are Jordan curves. Jonsson 
mentions they are in fact quasicircles, and shows, using Corollarv l2.4l that Cantor 
skews are Axiom A. Note such a map has only one basic set, which is topologically 
a Cantor set crossed with a circle. 

For this map, we established Axiom A with boxes from a (2 7 ) 2 x (2 5 ) 2 grid on 
23 = [— 2.1,2. 1] 2 x [— 1.28, 1.28] 2 . The box chain transitive component containing 
J j was of size (V; E) = (39,880; 1,447,640). We established box vertical expansion 
on this component by L = 1.125, with metric constants in [0.0024, 1], with average 
0.00529. See Figure Q] 

Example 5.3 (Cantor cross Basillica perturbation). One can generalize the above 
example by examining J p Cantor, and q z a small perturbation of any connected, 
hyperbolic polynomial, like the Basillica: w i— ► w 2 — 1. 

For the specific example f(z,w) — (z 2 — 6, w 2 — 1 — z/100), we achieved separation 
of the chain transitive components using boxes from a (2 10 ) 2 x (2 7 ) 2 grid on 23 = 
[-3.1, 3.1] 2 x [-1.732, 1.732] 2 (see Figure |2J|, but were unable to show Axiom A at 
that level, or even for a refinement of boxes from a (2 10 ) 2 x (2 8 ) 2 grid on 23. Trying 
to refine further maxed out our memory resources of 2GB of RAM. 

^All computations were performed on the Indiana University IBM Research SP system. Hence, 
this work was supported in part by Shared University Research grants from IBM, Inc. to Indiana 
University. 
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FIGURE 1. Let f(z,w) = {z 2 + 2, w 2 + z/10), as in Example I5~2l 
On the left is the box cover of J p , with boxes from a (2' ) 2 grid on 
[— 2.1,2.1] 2 . On the right is the slice of the box chain transitive 
components in the fiber over a repelling fixed point of p, z = 0.5 + 
1.32288z (necessarily in J p ), with boxes in the w-coordinate from 
a (2 5 ) 2 grid on [—1.28, 1.28] 2 . Since q z (w) is close to w 2 , there are 
two components: one containing the quasicircle J z , and the other 
a collection of boxes near the origin. 



Figure 2. For the map f(z, w) = (z 2 - 6, w 2 - 1 - z/100) of Ex- 
ample EH Jp is a Cantor set, while each q z is close to a Basillica. 
Shown above are slices of the box chain transitive components in 
two fibers: left over z — —3 = —(3, right over z = 3 = j3. Both show 
the same window, with boxes from a (2 7 ) 2 grid on [—1.732, 1.732] 2 
in the w-planc. Different box chain transitive components are 
colored different shades. There is one Basillica-like component 
containing Jj p , a second component near the attracting 2-cycle 
of the Basillica, {0,-1}, and a third component of islands skip- 
ping around the second component. This last component contains 
psewiio-rccurrent points only, and disappears for smaller boxes. 



Example 5.4 (Jonsson's nonproduct). The map f(z,w) — (z 2 — 90, w 2 + (z + 9)/4) 
is an example of the class of map which Jonsson showed to be Axiom A, but not 
conjugate to any product. J z in the fiber over the a-fixed point of p is a circle, 
while J z over the /3-fixed point of p is a Cantor set. 

For this specific map, we used Hypatia to prove Axiom A with boxes from a 
(2 7 ) 2 x (2 7 ) 2 grid on £ = [-10.1, 10.1] 2 x [-2.842, 2.842] 2 . Since J p is a Cantor 
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FIGURE 3. For Example 15.41 the map (z, w) — ► (z 2 — 90, w 2 + (z + 
9)/4). shown are box chain transitive components in four fibers: 
left to right over z = — 10 = — /3, z = —9 = a, z = 9 = —a, and 
z = 10 = /3, with boxes from a (2 7 ) 2 w-grid on [-2.842, 2. 842] 2 . 
The same window is shown for each fiber. 



set, there is only one box chain transitive component for which we must estab- 
lish vertical expansion: the one containing Jj . For a component of size (V; E) = 
(5,640; 209,572), we showed box vertical expansion by L = 1.25, with metric con- 
stants in [0.066, 1], with average 0.125. See Figure|31 

Example 5.5 (Generalization type 1 of Jonsson's nonproduct). The map f(z, w) = 
(z 2 — 90, w 2 + z/6 + 1.4 + 0.75i) satisfies Proposition 13.11 For this map, J z in 
the fiber over the a-fixed point is the Rabbit, while J z in the fiber over the f3- 
fixed point is a Cantor set. Note since A p is empty, we need only test vertical 
expansion on the box chain transitive component containing Jj p . For this map, 
we used Hypatia to prove Axiom A, with boxes from a (2 10 ) 2 x (2 9 ) 2 grid on CB = 
[-10.1, 10. 1] 2 x [-2.426, 2.426] 2 . 

See Figure 0] Note that the attracting 3-cycle of the Rabbit yields a saddle 3- 
cycle in C 2 , lying in the fiber over the a fixed point . This 3-cycle is contained in its 
own box chain transitive component , which must be separated from the component 
containing Jj before vertical expansion can be established. The component of Jj , 
of size (V; E) = (78,994; 2,066,558), is box vertically expansive by L = 1.0625, with 
metric constants in [0.000853, 1], with average 0.00315. 

Example 5.6 (Generalization type 2 of Jonsson's nonproduct). Example 13.21 is 
the simplest map of this type, with J p a real Cantor set, and the fibers containing 
(among other things) circles and Basillicas. We tested Hypatia on the map /(z, w) = 
(z 2 — 90, w 2 — z/20— 1/2) of this type. However, we were unable to achieve separation 
of all of the chain transitive components, hence could not test for Axiom A. The 
obstacle was memory usage. The deepest level we could completely build used 
boxes from a (2 11 ) 2 x (2 8 ) 2 grid on [-10.1, 10. 1] 2 x [-1.7203, 1.7203] 2 , but for these 
boxes the critical point of q z was still in the box chain transitive component for 
Jj p , and we used our entire 2GB available trying to build the next level. 

5.1. Conclusions and remaining questions. In conclusion, we would like to 
summarize our results, and provide the reader with questions for further study. 
We have begun to examine what kind of dynamics can occur for Axiom A skew 
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FIGURE 4. For Example 15. 51 the map (z,w) — ► (z 2 — 90, ■u; 2 -|-2:/6-|- 
1.4 + 0.75i), shown are box chain transitive components in four 
fibers: left to right over z = —10 = — (3, z = —9 = a, z = 9 = —a, 
and z = 10 = /?, with boxes from a (2 9 ) 2 w-grid on [-2.426, 2.426] 2 . 
The same window is shown for each fiber. 



products, by developing and implementing rigorous computer algorithms into a 
program Hypatia, to test whether a skew product of the form (z 2 + a, w 2 + bw + 
cz + e) is Axiom A. We have used Hypatia to verify the Axiom A property for 
Heinemann's and Jonsson's Examples; and further, derived formulas for two new 
classes of proposed Axiom A maps (generalizing Jonsson's examples), and used 
Hypatia to prove Axiom A for some examples in our first new class. 

Polynomial skew products are a rich area for study. Wc would like to encour- 
age the interested reader to begin an exploration of this family, by providing the 
following list of questions which arose during our recent introduction to this area. 

Problem 5.7. Produce a combinatorial description of the dynamics of maps sat- 
isfying Proposition or of the same type as Examples 13.21 and \'6.'6\ For example, 
by employing the theory of "random" dynamical systems, and studying a random 
composition of two hyperbolic polynomials of C (like z 2 and z 2 — 1 for Examplc l3.2(l . 

Question 5.8. If a polynomial skew product is an Axiom A stable perturbation 
of a product, then arc all of the fibers homcomorphic? 

Question 5.9. Is there an Axiom A polynomial skew product with J p connected 
but whose fiber Julia sets J z for z £ J p are not all homcomorphic? Jonsson (\12\. 
Example 9.6) and Proposition 13 . II both provide examples of Axiom A maps whose 
fibers are not homcomorphic, but J p is a Cantor set. 

Question 5.10. What can be said about polynomial skew products using the 
techniques of holomorphic motions, in dynamical space and/or parameter space? 

Question 5.11. What further measure theoretic results hold for skew products? 
For example, is there a "typical" J Zl for z € J pi in the measure theoretic sense? 

Question 5.12. How can the theory of partial hypcrbolicity, which is well de- 
veloped in the real variables setting, and can pertain to skew products of real 
manifolds, shed light on the structure of polynomial skew products of C 2 ? 
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Question 5.13. Jonsson (|12j. Theorem A. 6) establishes structural stability on 
Jj p , for Axiom A polynomial skew products. When (if ever) can structural stability 
be extended to hold over a larger set? 

Appendix A. Rigorous Arithmetic 

On a computer, we cannot work with real numbers; instead we work over the 
finite space F of numbers representable by binary floating point numbers no longer 
than a certain length. For example, since the number 0.1 is not a dyadic rational, 
it has an infinite binary expansion. Thus the computer cannot encode 0.1 exactly. 
Interval arithmetic (I A) provides a method for maintaining rigor in computations, 
and also is natural and efficient for manipulating boxes. The basic objects of I A 
are closed intervals, [a] = [a, a] G IK, with end points in some fixed field, K. An 
arithmetical operation on two intervals produces a resulting interval which contains 
the real answer. For example, 

[a] + [b] := [a + b,a + b\ and [a] - [b] := [a - b, a - 6] 

Multiplication and division can also be defined in IA. 

Since an arithmetical operation on two computer numbers in F may not have a 
result in F, in order to implement rigorous I A we must round outward the result of 
any interval arithmetic operation, e.g. for [a], [b] G IF, 

[a] + [b] := [ia + M ,t« + &T] . 

where [xl denotes the largest number in F that is strictly less than x (i.e., x 
rounded down), and denotes the smallest number in F that is strictly greater 
than x (i.e., x rounded up). This is called IA with directed rounding. 

For any let Hull(x) be the smallest interval in F which contains x. That is, 

if x G F, then Hull(a;) denotes [x,x\. If x G R \ F, then Hull(x) denotes [[x[ , ]x]]. 
Similarly, for a set S C R, we say Hull(iS') for the smallest interval containg S. 
Whether Hull(5) is in IR or IF should be clear from context. 

In higher dimensions, IA operations can be carried out component-wise, on in- 
terval vectors. So if x G R™, then Hull(x) = Hull(xi) x ••• x Huil(z n ), and if 
S C R™, then Hull (5) is the smallest vector in IF" (or IR") containing S. Note to 
deal with intervals in C™ we simply identify C™ with R 2n . Thus a box in C = R 2 
is an interval vector of length two. Our extensive use of boxes is designed to make 
I A calculations natural. 

To compute the image of a point x under a map / using IA, first convert x to 
an interval vector X = Hull(a;), then use a (carefully chosen) combination of the 
basic arithmetical operations to compute an interval vector F(X), such that we are 
guaranteed that Hull(/(x)) C F(X). If /: R™ —> R™ is continuous, then an interval 
extension of /, F = F(f), is a function which maps a box B in W L to a box F(B) 
containing f(B), i.e., F(B) D Hull(/(B)). Usually, we would like F(B) to be as 
close as possible to Hull(/(-B)). We shall not discuss here how to find the best F. 

Each time an arithmetical calculation is performed, one must think carefully 
about how to use IA. For example, IA is not distributive. Also, it can easily create 
large error propagation. For example, iterating a polynomial map on an interval 
vector (which is not very close to an attracting period cycle) will produce a very 
large interval vector after only a few iterates. That is, if B is a box in R 2 = C, and 
one attempts to compute a box containing p*°(B), for p c (z) = z 2 + c, by: 
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for j from 1 to 10 do 
' B= Pc (B) 

then the box B will likely grow so large that its defining bounds become machine 
oo, i.e., the largest floating point in F. 

For further background on interval arithmetic, see We use IA for all 

the rigorous computations in the computer program Hypatia. The IA routines were 
all provided by the PROFIL/BIAS package, available at [T%] . 
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